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We study how the Fermi energy of a graphene monolayer separated from a conducting substrate 
by a dielectric spacer depends on the properties of the substrate and on an applied voltage. An 
analytical model is developed that describes the Fermi level shift as a function of the gate voltage, of 
the substrate work function, and of the type and thickness of the dielectric spacer. The parameters of 
this model, that should describe the effect of gate electrodes in field-effect devices, can be obtained 
from density functional theory (DFT) calculations on single layers or interfaces. The doping of 
graphene in metal | dielectric | graphene structures is found to be determined not only by the difference 
in work function between the metal and graphene and the dielectric properties of the spacer but 
potential steps that result from details of the microscopic bonding at the interfaces also play an 
important role. The doping levels predicted by the model agree very well with the results obtained 
from first-principles DFT calculations on metal | dielectric | graphene structures with the metals Al, 
Co, Ni, Cu, Pd, Ag, Pt or Au, and a /i-BN or vacuum dielectric spacer. 

PACS numbers: 73.63.-b, 73.40.-c, 73.22.Pr, 73.20.Hb, 81.05.U- 



I. INTRODUCTION 

The successful extraction of graphene, single monolay- 
ers of graphite, has led to the discovery of new phys- 
ical phenomena that result from the unique electronic 
structure of this two-dimensional system.^ Many of the 
remarkable properties of graphene have been explored 
using field effect devices^'^ in which a graphene flake is 
placed on top of an oxidized silicon wafer and contacted 
with metal (source and drain) electrodes. 

Direct contact between graphene and a metal electrode 
generally results in a transfer of electrons between the 
two materials. Depending upon the metal species used, 
the contact causes graphene to be doped n- or p-type in 
equilibrium.^'^ Metal-contact doping also induces a lat- 
eral space charge and corresponding band bending in the 
graphene adjacent to, but not covered by, the metal con- 
tact. The two-dimensional character of graphene and the 
low density of states (Dos) around its Fermi level then 
result in an extended space charge region.^ The space 
charge distribution is inhomogeneous, which can lead to 
p-p', n-n', p-n or, with gating, to even more complex 
profiles such as p-n-p' within the space charge region. 

In this paper we focus on electrostatic doping, where 
application of a voltage difference between graphene and 
a back gate changes the position of the Fermi energy 
with respect to the Dirac point^'^"^^ and a nonlinear 
shift of the graphene Fermi energy as a function of the 
gate voltage has been observed in scanning tunneling 
spectroscopy (STS) experiments^^ and in work function 
measurements.^^ In experiments such as these, where a 
Si02 dielectric was used with a doped Si gate, it is be- 
lieved that intrinsic inhomogeneities in Si02 as well as 
surface roughness are reasons that much lower charge 
carrier mobilities are observed than for freely suspended 



graphene. 

In this respect, a much better gate dielectric for 
graphene based electronics is hexagonal boron-nitride, 
/i-BN.^^'^^^^ Graphene adsorbed on /i-BN has been 
shown^^^^ to be much flatter than on Si02 and to have 
highly mobile charge carriers. This favorable situation 
can be related to the fact that /i-BN has a layered crystal 
structure similar to that of graphite; it consists of atoms 
that are strongly bonded in plane in a honeycomb lat- 
tice, whereas the binding between these planes is much 
weaker. But where graphene is a semi-metal, /i-BN is 
an insulator. Graphene binds weakly to /i-BN with an 
interaction similar to that between /i-BN monolayers, or 
between the graphene layers in graphite. All of this 
means that adsorbing graphene onto /i-BN perturbs its 
unique electronic structure minimally. 

Like (multilayer) graphene, layers of /i-BN can be pre- 
pared by micromechanical cleavage, and can be trans- 
ferred to a metal substrate to form a gate electrode. ^^'^^ 
Alternatively, single layers of /i-BN have been grown di- 
rectly on the (111) surface of transition metals such as 
Cu, Pt, Rh, Pd and Ni.^^ ^e Graphene |/i-BN heterostruc- 
tures have been grown on Ni(lll)^^ and more recently on 
Ru(OOOl).^^ Very recently, structures of alternating /i-BN 
and graphene layers have been constructed and shown to 
operate as a field-effect tunneling transistor. 

Obviously the properties of the gate electrode play an 
important role in determining the behavior of graphene 
field effect devices. In this paper we develop an analytical 
model for the electrostatic doping of graphene as a func- 
tion of the electric field generated by a gate bias and test 
the model against results obtained using first-principles 
density functional theory (DFT) calculations on metal 
gate I dielectric I graphene (M|D|Gr) stacks in the presence 
of an external electric field across the stack. The results 
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of calculations for the particular case of M = Cu were 
recently briefly reported.^'' 

The electrostatic model for the M|D|Gr structure is 
based upon a parallel plate capacitor geometry. The in- 
put parameters to the model are the work function of 
the metal gate, the dielectric constant and thickness of 
the gate dielectric, and the potential steps formed at the 
interfaces between the metal and the dielectric, and be- 
tween the dielectric and the graphene sheet. Such poten- 
tial steps originate from interactions between the two ma- 
terials at an interface, and they can make a no n- negligible 
contribution to the potential profile. ^'^'^^'^^'^^ All the pa- 
rameters are accessible to DFT calculations on single lay- 
ers or interfaces. 

In the DFT calculations on the full M|D|Gr structures, 
we use primarily /i-BN as dielectric material and study 
the series of metals Al, Co, Ni, Cu, Pd, Ag, Pt, and Au, 
whose work functions span a considerable range, from 4.2 
to 6.0 eV. /i-BN is physisorbed on Al, Cu, Ag, Pt, and Au 
and chemisorbed on Co, Ni, and Pd allowing us to study 
the effect of the bonding between metal and dielectric 
on the electrostatic doping of graphene. To assess the 
influence of the dielectric, we examine the effect of using 
a vacuum "dielectric" instead of /i-BN. 

The paper is organized as follows. In the next section 
we formulate the analytical model for electrostatic dop- 
ing of graphene. We describe the DFT calculations and 
the (electronic) structures of the metal |/i-BN| graphene 
stacks in Sec. III. The results of the calculations are pre- 
sented in Sec. IV and compared to the predictions of the 
model. Finally, Sec. V presents a summary plus some 
conclusions. 



II. ELECTROSTATIC DOPING MODEL 

We model the M|D|Gr structure as a parallel plate ca- 
pacitor. Such a model has been used successfully to ex- 
plain the Fermi level shift in graphene directly absorbed 
on a metal substrate. ^'^ Inserting a dielectric is expected 
to modify the potential profile as shown schematically in 
Fig. 1. The shift of the Fermi level, AEp^ with respect to 
the charge neutrality point in graphene, i.e. the conical 
points of the band structure, is given by 



Evac*-7<: 



AEf = W-WGr 



(1) 



with Wcr and W the work functions of graphene in vac- 
uum, and of the complete M|D|Gr stack, respectively. 
Note that a negative value of A£^f represents electron 
doping, a positive value hole doping. W can be written 



W = Wm- AF, 



(2) 



with Wm the work function of the clean metal substrate, 
and AV the potential difference across the stack. The 
latter will in general be a function of the dielectric thick- 
ness d. 
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FIG. 1. (color online) Schematic drawing of the 
metal I dielectric I graphene (M|D|Gr) structure. The electro- 
static potential profile in the structure is indicated by the 
wide dashed (orange) line at the top. AV is the shift of the 
work function of the structure with respect to that of the 
clean metal, VKm; Am|D5 A^iGr are the potentials steps at the 
interfaces; A^f is the Fermi level shift in graphene caused by 
the electrostatic doping. 



We model the potential difference as 
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where Ei^t is the electric field inside the dielectric, — e 
is the charge of an electron, and Am|d and A^iGr 
are potential steps at the metal | dielectric and the 
dielectric I graphene interfaces, respectively. Interface 
bonding between two materials usually results in the for- 
mation of a dipole layer that is accompanied by a poten- 
tial step^^^^ localized at the interface. Interface dipoles 
are intrinsic properties of interfaces and even if the in- 
teraction at an interface is as weak as in the case of ph- 
ysisorption, they can be sizable. We will assume that 
interface potential steps are independent of applied elec- 
tric fields for field magnitudes of practical interest and 
test this assumption with DFT calculations. 

In the absence of a dielectric spacer, an externally ap- 
plied gate voltage gives rise to an electric field £^ext- In 
the following, it will be convenient to use £^ext as an inde- 
pendent variable. The internal field, Eq. (3), is obtained 
from the standard electrostatic boundary condition 



eo (^ext - I^Eint) = Cr, 



(4) 



with hz the relative dielectric constant of the dielectric 
layer, and a the surface charge density of the graphene 
sheet. 

The latter can be related to the Fermi level shift AE-p 
by noting that charge in graphene is introduced by oc- 
cupying states starting from the charge neutrality point, 
a = e J^^^ D{E)dE, where the charge neutrality point 
corresponds to the conical points in the graphene band 
structure. The density of states near the conical points is 
described well by the linear function, D{E) = \E\ Dq/A, 
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with Do = 0.102/(eV^ unit cell), and A = 5.18 as the 
area of a graphene unit cell, so that 



<T = sign(A^F)^(AEF)2. 



(5) 



The value of I^o is sensitive to the k-point grid used 
for the Brillouin zone sampling, and it depends on the 
substrate used; technical details are discussed in the 
Appendix.^^ 

Equations (l)-(5) can be combined to give a closed 
expression for AE-p in terms of the applied external field 

^ext 



A^F(^ext) = ± 



+ 2Doa{d/K)^ |e(^ext + ^o)| - 1 



Doad/hi 



(6) 



where a = e'^ / {cqA) = 34.93 eV/A for a 1 x 1 graphene 
unit cell. The sign of AE-p is determined by the sign of 
Ee^t +£^0, where 



ae 



with 



Vo = Wm- Wcr - A 
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(7) 



(8) 



Eq is the value of the intrinsic electric field across the 
dielectric in the absence of any external electric field, see 
Fig. 1, and it can be substantial when d is sufficiently 
small, i.e., for thin dielectric spacers. It results from an 
electron transfer in equilibrium between the metal and 
graphene across the dielectric spacer that leads to an 
intrinsic doping of graphene. A compensating external 
field £^ext = —^0 is required to bring the Fermi level to 
the conical points and make graphene electrically neutral, 
AEf = 0. 

According to Eqs. (l)-(3), the Fermi level shift AE-p 
and the internal electric field Ei^t are linked by a linear 
relation. This means that both these quantities are simi- 
lar non-linear functions of the external electric field £^ext , 
Eq. (6), whose non-linearity is determined by the doping 
charge on graphene. Figure 2 shows ^int as a function of 
E'ext for two different values of the dielectric constant hz. 
The points in these curves where the derivative is maxi- 
mal correspond to undoped graphene, i.e., AE-p = and 
Eext = —^0- It follows from Eq. (4) that at these points 
the curves cross the lines Eq^^/k^ see Fig. 2. From Eq. 
(7) it follows that the corresponding internal electric field 
is Eint = -V^/de. 

The £^int curves cross at £^int = 0. From Eqs. (l)-(5) 
it is easily seen that this occurs when ^ext = ^io, where 
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FIG. 2. (color online) The internal field ^int (solid lines) 
inside the dielectric spacer of Fig. 1, as a function of the 
applied external field ^ext, for two values of n. The black line 
corresponds to a Cu|/i-BN| Graphene structure with 5 layers 
of /i-BN and n — 2.8. The red line would be the response of 
the same system, but with k = 1.0 instead of 2.8. The curves 
cross at E^^t — ^io, see Eq. (9), where ^int = 0. At fields 
^ext = —^0, see Eq. (7), the curves cross the lines Ee^t/i^ 
(dashed lines). 



is independent of d and z^, i.e., all curves A£^F(^ext) for 
different dielectric thicknesses cross at AE^ = Vq for 
£^ext = ^io- This crossing point could be used to deter- 
mine Vo from experiment. 

Equation (6) constitutes the model we will use in the 
remainder of this paper. Assuming that the properties of 
graphene are known, the model depends on the parame- 
ters Wm^ ^m|D7 ^D|Gr7 and k,. All of these parameters 
can be obtained from DFT calculations. The thickness 
of the dielectric spacer d is somewhat ill-defined at the 
atomic level. As in Ref. 5, we set d = D — do^ where D 
is the distance from the center of the top layer of metal 
atoms to the graphene plane, and do = 2.4 A. 

The model can also be expressed in terms of an applied 
gate voltage, writing the latter as 



Vg = -eEe^td/K.. 
Equation (6) then becomes 



(11) 



^^^^^VI±^^D^MEM^, (12) 

aDo d/ Hi 

and the sign of AE^ is determined by the sign ofVo — Vg. 

III. COMPUTATIONAL DETAILS 

In this section we briefiy describe the most salient fea- 
tures of the first-principles DFT calculations that we use 
to determine the parameters entering the model defined 
in the previous section. We calculate electronic ground 
state energies within the framework of the plane-wave 
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FIG. 3. (color online) Top views of the energetically favorable 
configurations of /i-BN on a) Cu, Ni(lll), Co(OOOl), and b) 
Pt, Pd, Au, Ag, Al (111) in 1x1 and ^3x^3 surface unit 
cells, respectively. 



PAW pseudopotential method, at the level of the lo- 
cal density approximation (LDA),^^ as implemented in 
VASP^^^^ and use them to optimize the geometries of 
M|D|Gr structures for a number of representative metals. 
The structures are modeled as periodic slabs separated 
by a vacuum region ~ 15 A thick. A dipole correction is 
applied to avoid spurious interactions between periodic 
images of the slab.^^ A plane wave kinetic energy cutoff 
is set at 400 eV. 

A 36 X 36 k-point grid is used to sample the surface 
Brillouin Zone (BZ) and structural relaxations are per- 
formed using a smearing method for the BZ integration,^^ 
with a smearing parameter of 0.2 eV. Total energies and 
charge densities of optimized structures are then recalcu- 
lated using the tetrahedron scheme. The special points F, 
K' and M of the graphene band structure are explic- 
itly included in the BZ sampling, which is important for 
an accurate description of the electronic structure around 
these points. 

The electronic self-consistency criterion is set to 10~^ 
eV. Such a strict criterion is required to assure conver- 
gence of the charge distribution. We consider only close- 
packed metal surfaces with three-fold rotational symme- 
try, i.e., the (111) surface for fee metals, and the (0001) 
surface for hep metals. A slab of six layers of atoms is 
used to represent the metal surface. The /i-BN layers are 
adsorbed on one side of the slab to represent the dielec- 
tric, and graphene is adsorbed on top of this. To model 
the effect of the gate potential, a sawtooth-shaped poten- 
tial is applied across the whole slab.^^ The carbon atoms 
in graphene, the boron and nitrogen atoms in the /i-BN 
slab, and the separation between the top two atomic lay- 
ers of the metal were allowed to relax during the geome- 
try optimization until the total energies were converged 
to within lO"'^ eV. 

The LDA is found to give a reasonable description of 
the weak interaction between /i-BN layers and between a 
/i-BN layer and graphene. Likewise, LDA describes the 
interaction between graphene and metal surfaces reason- 
ably well.^^ We expect a comparably good description 
of the interaction between /i-BN and metal surfaces and 
use the LDA functional throughout. 




FIG. 4. (color online) Side views of the energetically favorable 
configurations of metal |/i-BN| graphene structures, a) Cu, Ni, 
Co and b) Au, Pt, Ag, Pd, AL 



For metal|/i-BN interface calculations, we fixed the lat- 
tice constant of /i-BN to its optimized LDA value of 
<^hex = 2.49 A. The in-plane lattice parameter of the 
metal surfaces was then adjusted so that for Cu, Ni and 
Co we could use a common 1x1 surface unit cell. Fig. 3a. 
For Au, Pt, Al, Ag, and Pd, a a/3 x a/3 surface unit cell 
matches very well to a 2x2 surface unit cell of /i-BN, 
Fig. 3b. Adopting this to model the interface involves 
changing the in-plane lattice parameter of the metals by 
less than 2% compared to the optimized values. Such a 
change does not greatly perturb the electronic structure 
of the metals, especially after allowing the top two metal 
layers to relax in a direction perpendicular to the surface. 
The small change in the metal lattice parameter is found 
to have only a small effect on the properties we are inter- 
ested in. For instance, the work function of the Cu(lll) 
surface is increased by a mere 0.05 eV and the density of 
states at the Fermi level is unaltered. 

In the most stable configuration of a /i-BN sheet on Cu, 
Ni and Co, nitrogen atoms are adsorbed on top of metal 
surface atoms and boron atoms are at hollow sites as 
shown schematically in Fig. 3a. The most stable struc- 
ture for /i-BN on Au, Pt, Al, Ag, and Pd is shown in 
Fig. 3b. In an interface unit cell containing 3 surface 
metal atoms, 4 boron and 4 nitrogen atoms, one boron 
and one nitrogen atom is absorbed on top of a metal 
surface atom; all other boron and nitrogen atoms are on 
bridge sites. These structures are used to calculate the 
metal|/i-BN interface dipoles. 

For calculations on metal |/i-BN| graphene stacks we use 
the structures shown schematically in Fig. 4. Here we 
choose the lattice constant of graphene equal to its op- 
timized LDA value a = 2.445 A, and adapt the in-plane 



5 




r K Mr KM 



FIG. 5. (color online) Left panel: the band structure of a 
Au(lll)|/i-BN monolayer I graphene stack. The eigenstates of 
graphene and /i-BN with strong tt character are highlighted in 
black and green (gray), respectively where the symbol sizes 
represent the amount of pz character. The plotting order 
is total, graphene pz, /i-BN pz. Right panel: the superim- 
posed band structures of free-standing graphene and a /i-BN 
monolayer. The graphene (/i-BN) bands have been shifted up 
(down) by 0.08 (1.28) eV. 



lattice constants of /i-BN and the metal accordingly. This 
means that the /i-BN is compressed in plane by 1.8% and 
the metal by up to 3.8%. The same procedure was used in 
our earlier studies of graphene on metal surfaces^ and on 
bulk /i-BN,^^ the basic rationale being that a change of 
the in-plane lattice parameter of the metal or of /i-BN will 
have a smaller effect on their electronic structures (at the 
Fermi energy) than a corresponding change in graphene 
would have on its properties. For instance, compressing 
the Cu(lll) lattice by 4% increases its work function by 
only 0.08 eV. As a second example, the potential step 
at the Cu(lll)|/i-BN interface is Acu|bn =1.12 eV for 
<^hex = 2. 445 A, the lattice constant of graphene, and 1.18 
eV for tthex = 2.49A, the h-BN lattice constant. 

/i-BN layers are stacked in AA' fashion, with boron 
on top of nitrogen and vice versa, see Fig. 4a. The 
graphene layer is placed on top of /i-BN as in Ref. 15. 
The optimized LDA interlayer separation in the /i-BN 
slab is 3.24 A, which is similar to what is found in bulk 
/i-BN. For the Au(lll)| monolayer /i-BN| graphene struc- 
ture (Fig. 4b) the calculated equilibrium separation be- 
tween the top /i-BN layer and graphene is 3.21 A, which 
is similar to the calculated separation of 3.22 A between 
graphene and bulk /i-BN,^^ as expected. 

In the analytical model outlined in Sec. II we have as- 
sumed that, apart from a shift in the Fermi level, the elec- 



tronic structure of graphene in a met al|/i-BN| graphene 
stack is unchanged compared to a free-standing sheet of 
graphene. We can check this statement explicitly by ana- 
lyzing the electronic structure obtained from a DFT cal- 
culation on a M|/i-BN|Gr stack. As an example, the left 
panel of Fig. 5 shows the band structure of a Au(lll) sub- 
strate with a monolayer of /i-BN and graphene on top. 
The graphene and /i-BN pz character of the eigenstates is 
indicated to help identify the bands. For comparison, the 
electronic structures of free-standing graphene and of a h- 
BN monolayer are shown superimposed in the right panel 
of Fig. 5. To make the Dirac points of graphene coincide 
in energy, the band structure of free standing graphene 
in the right panel has been shifted downwards by 0.08 
eV corresponding to the doping level of graphene in the 
Au|/i-BN monolayer I graphene structure. The top of the 
/i-BN valence band, at K, is positioned at —1.28 eV, cor- 
responding to the Schottky barrier for holes found for the 
Au|/i-BN structure. Apart from these energy shifts, the 
band structures of both graphene and /i-BN in the Au\h- 
BN I graphene stack are hardly changed with respect to 
their free-standing counterparts. 

This reflects the weakness of the interaction between 
Au and /i-BN, and between /i-BN and graphene. Some 
metal substrates have a stronger interaction with /i-BN 
which does affect the /i-BN derived bands. However, as 
far as the doping levels of graphene are concerned, the 
effects of that interaction can be fully incorporated in the 
potential step Am|bn associated with the metal |/i-BN in- 
terface that will be discussed in the next section. If the 
/i-BN and the graphene lattices are assumed to be com- 
mensurate, the sublattice-symmetry-breaking potential 
of /i-BN induces a small band gap of ~ 54 meV between 
the tips of the Dirac cones; apart from the formation of 
this small gap, the DoS remains linear in energy. If the 
lattices are incommensurate, this gap disappears. ^^^^'^^ 
In the following this point will not be important. 



IV. RESULTS 

Conceptually, the simplest model we can consider for 
a study of electrostatic doping of graphene based upon 
first-principles calculations is a metal | vacuum | graphene 
stack. Fig. 6 shows the shift of the Fermi level in 
graphene calculated using Eq. (1), as a function of the ex- 
ternal electric field for the Cu(lll) | vacuum |graphene het- 
erostructure, where the vacuum "thickness" d = 7.18 A. 
The points in this figure are the calculated DFT values, 
the solid line represents the model described by Eqs. (6) 
and (7) with Hi = 1 and Am|d = ^DiGr = for the 
vacuum dielectric. The match between the model (solid 
curve) and the DFT values (squares) is excellent. 

At zero external bias, the work function difference be- 
tween the Cu(lll) surface and graphene sheet transfers 
electrons from one to the other until equilibrium is es- 
tablished. The transferred charge sets up an intrinsic 
electric field £^0 that is determined by the work function 
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FIG. 6. (color online) The Fermi level shift versus external 
electric field strength for a Cu(lll)|/i-BN|graphene structure 
with two layers of /i-BN (black) and for the same structure 
with the h-BN layers replaced by vacuum (red). The lines 
represent the model described by Eq.(6) with k = 2.80 and 
1.0, respectively. The symbols are the values calculated with 
DFT. 



the range 0.05-0.45 V/A, we find k = 2.80 ± 0.02. In our 
modelling in the remainder of this paper, we use Hi = 2.80. 
Fig. 6 shows that with this value of the shape of the 
curve gives an excellent description of the DFT results 
for the Cu(lll)|/i-BN|graphene heterostructure. 

The shift between the two curves in Fig. 6 results from 
potential steps Am|d and A^iGr that are formed at the 
metal|/i-BN and /i-BN|graphene interfaces, respectively. 
These potential steps are sufficiently large not only to 
compensate for the difference in work functions of the 
Cu(lll) surface and free-standing graphene, but even to 
reverse the intrinsic electric field (Eq. (7)). It means 
that in the Cu(lll)|/i-BN|graphene stack at zero bias, 
graphene is n-type doped, and one needs to apply a pos- 
itive external field to undo the graphene doping and en- 
force AE-p = 0. Because the interface potential steps play 
such an important role, we study them in more detail in 
the next subsection. 



Potential steps at metal |/i-BN and 
^-BN I graphene interfaces 



difference, see Eq. (7). The work functions of the clean 
Cu(lll) surface and free-standing graphene calculated in 
the LDA are 5.25 eV and 4.48 eV, respectively, which in- 
dicates that at zero external bias electrons are transferred 
from graphene to Cu so that graphene is doped p-type. 
A compensating negative external field can be applied 
to undo the graphene doping so AEip = 0. The values 
given by the model of the compensating field and the 
Fermi level shift at zero bias are in excellent agreement 
with the DFT results on the Cu(lll) | vacuum |graphene 
system, see Fig. 6. 

Insertion of a /i-BN dielectric spacer between the Cu 
and graphene yields a system that can be realized in prac- 
tice. The Fermi level shift in graphene for a Cu(lll)|/i- 
BN I graphene structure with two layers of /i-BN (corre- 
sponding, by construction, precisely to the vacuum thick- 
ness used in the previous calculation) is also shown in 
Fig. 6 as a function of the external electric field. Com- 
paring the results obtained with a /i-BN spacer and a 
vacuum spacer, two qualitative differences are immedi- 
ately obvious. The curve calculated with /i-BN results is 
shifted to the right, and it is "flatter". 

The latter effect results largely from the dielectric con- 
stant of the spacer; Eq. (6) shows that AE^ decreases 
with increasing dielectric constant. We determined the 
dielectric constant of /i-BN from calculations on very thin 
/i-BN films as the ratio between an applied external elec- 
tric field and the internal electric field inside the slab"^^ 
determined from the shift of the Is Nitrogen core levels. 
Because these levels are localized on each layer, their en- 
ergies shift according to the local electrostatic potential. 
The shift is linear in the position of the layer, and the 
slope yields the internal electric field. From calculations 
on /i-BN slabs consisting of 2-10 layers in electric fields in 



In general, potential steps are formed at interfaces be- 
tween different materials. They result from interface 
dipole layers, caused by chemical interactions and/or 
charge transfer between two materials comprising an in- 
terface. Substantial interface dipoles can be formed even 
if the interaction between the two materials is relatively 
weak as happens when noble gas atoms are adsorbed 
on, or closed-shell molecules are physisorbed on metal 
surfaces. ^^'^^'^^'^^ Even in the absence of a direct inter- 
action, the difference between the work functions of two 
metals leads to charge transfer and the formation of con- 
tact potentials. 

The potential step Am|d at the metal | dielectric inter- 
face can be obtained by comparing the work function Wu 
of the clean metal surface without the adsorbate to that 
of the metal with the dielectric adsorbed on it, I^m|d- 
Similarly, the potential step A^iGr is given by the dif- 
ference between the work functions on the dielectric side 
and the graphene side of a D|Gr slab 

Am|D =Wm- Wm|d; AoiGr = V^DIGr " V^Gr|D. (13) 

It turns out to be sufficient to study graphene adsorbed 
on two layers of /i-BN to find a converged value for 
Abnigf = 0-14 eV. Similarly, a single /i-BN layer ad- 
sorbed on a metal substrate is sufficient to determine 
Am|bn7 demonstrating the true interface nature of these 
potential steps. The values of Am|bn calculated for the 
metals considered in this paper, given in Table I, are 
seen to be quite substantial, ranging from 0.4 eV for 
Al(lll)|/i-BN to 1.8 eV for Co(0001)|/i-BN. In all cases 
Am|bn is positive, meaning that adsorption of /i-BN on 
a metal substrate lowers the work function. In other 
words, the interaction between /i-BN and the metal sub- 
strate leads to a net displacement of the electronic charge 
distribution into the metal. 
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TABLE L Calculated and experimental potential steps Am|bn 
at the metal|/i-BN interface, deq is the equilibrium separation 
of the /i-BN layer from the metal surface. The /i-BN lattice 
constant ahex = 2.49 A is used for the surface unit cells, ex- 
cept in the last two rows, where the graphene lattice constant 



Q^hex 


= 2.445 A is used. 






M 




eq \-t^) 




Co 


1.80 


1.92 


5.52 


Ni 


1.73 1.5-1.8"'^ 


1.96 


5.52 


Pd 


1.25 1.3 


2.47 


5.53 


Cu 


1.18 0.24^0.8-1.1^ 


2.97 


5.17 


Pt 


1.04 0.9' 


3.04 


5.98 


Ag 


0.83 


3.20 


4.83 


Au 


0.79 


3.24 


5.55 


Al 


0.41 


3.55 


4.25 


Or 


0.11 


3.21 


4.63 


Gr^ 


0.14 


3.22 


4.48 


Cu^ 


1.12 0.24^0.8-1.1^ 


3.11 


5.25 




"Refs. 45-47, 'Ref. 45, ^Ref. 25, 


^Ref. 44. 





The metal|/i-BN equilibrium distance deq indicates 
that there are two different adsorption regimes. Ad- 
sorption of /i-BN on the Co(OOOl), Ni(lll) and Pd(lll) 
surfaces leads to equilibrium distances d^^ < 2.5 A, val- 
ues which are typically associated with chemisorption, 
whereas the values of deq > 3.0 A for adsorption on the 
Cu, Pt, Ag, Au, and Al(lll) surfaces are typical of ph- 
ysisorption. 

The calculated values of Am|bn increase with decreas- 
ing separation of the /i-BN plane from the metal sub- 
strate, as shown in Table L A smaller separation leads 
to an increased metal |/i-BN interaction, and hence to a 
larger deformation of the charge density and a larger po- 
tential step at the interface. The calculated Am|bn gen- 
erally agree well with the experimental values. However, 
for Cu|BN, different experiments report quite different 
results. ^^'^"^ The values we calculate are in better agree- 
ment with values recently reported in Ref. 44. In view 
of the results shown in Table I, the experimental value 
^Cu|BN = 0-24 eV reported in Ref. 25 is strikingly low. 

An interface potential step can have different origins. 
Chemical interactions between two materials can result 
in an ordered interface layer of bond dipoles.^^'^^ Even 
when the interaction is weak, as in the case of physisorp- 
tion, it is argued that the Pauli exchange interaction be- 
tween the adsorbate and the substrate can distort the 
charge density distribution significantly and drive a net 
displacement of charge towards the "softer" material, in 
this case the metal substrate, which lowers the work func- 
tion. This is commonly called the "pushback" or "pillow" 
effect.43 

In addition to a potential step induced by a direct in- 
teraction at the interface, one can have a contribution 
originating from an electron transfer to equilibrate the 
Fermi level between the two materials on either side of 
the interface. Such a transfer obviously occurs only if 
the Fermi level is in a range where the density of states 



of both materials is non-zero. At interfaces between a 
metal and an organic molecular material, where charge 
is exchanged between strong donor or acceptor molecules 
and the metal surface, this leads to pinning of the Fermi 
level by the molecules. 

One does not expect this to be the case for /i-BN ad- 
sorbed on a metal surface. /i-BN is a large band gap 
insulator, and in most cases the Fermi level of the metal 
is well within the /i-BN band gap. The potential steps 
^M|BN and Abnigf can therefore be ascribed to the di- 
rect interaction between the materials at the interface. 

We can visualize the electronic displacement and dipole 
formation at the interface by looking at the electron den- 
sity of the entire system minus the electron densities 
of the two separate constituent materials. As only the 
component perpendicular to the interface is relevant in 
the present context, it is convenient to work with plane- 
averaged electron densities n{z) = J n{x,y^ z)dxdy /A in 
terms of which the electron charge density p = —en. The 
electron displacement in a metal |/i-BN structure is de- 
fined as 

AnM|BN(^,^ext) = 

^M|Bn(^, £^ext) - ^m(^, ^ext) " ^Bn(^, ^^ext), (14) 

where nM|BN is the electron density of the composite 
system, nu and ubn are the electron densities of the 
clean metal substrate and the isolated /i-BN layer, re- 
spectively, and we introduce an explicit £^ext dependence 
because we wish to study the dependence of the inter- 
face electronic displacement on the external electric field. 
Any electronic displacement An{z) can be related to a 
potential (energy) step A by simple electrostatics, i.e., 
A = -e^ / zAn{z)dz/eo. 

The calculated electronic displacement without an ex- 
ternal field, AnM|BN(^7 ^ext = 0)? is shown in Fig. 7 for 
four different metal |/i-BN interfaces. The amplitude of 
AnM|BN(^) is seen to decrease with increasing separa- 
tion of the adsorbed /i-BN layer. The oscillation pat- 
terns illustrate the two different kinds of behavior dis- 
cussed above. For /i-BN adsorbed on the Cu(lll) and 
Au(lll) surfaces, the pattern of AnM|BN(^) in the re- 
gion between the metal surface and the /i-BN layer is 
that of a simple dipole. It corresponds to an accumula- 
tion of electrons close to the metal surface and a deple- 
tion close to the /i-BN layer, consistent with the pillow 
effect discussed above. The simple interface dipole and 
the large metal |/i-BN equilibrium separations (Table I) 
found for these systems are in line with physisorption of 
/i-BN on these metals. In contrast, the oscillation pat- 
tern of AnM|BN(^) for /i-BN adsorbed on Ni(lll) and 
Pd(lll) surfaces is much more complicated. It is ac- 
companied by a small metal-/i-BN equilibrium distance, 
which is consistent with chemisorption. 

In the model described by Eqs. (6)- (8), it was as- 
sumed that the interface potential steps Am|d and A^iGr 
are independent of an applied external field. That im- 
plies that the corresponding interface dipoles should 
be independent of the external field. Figure 8 shows 
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FIG. 7. (color online) Plane- averaged electron density dif- 
ference AnM|BN(2;, 0) for Ni, Pd, Cu and Au (111) surfaces 
covered with a monolayer of /i-BN. The vertical dashed lines 
represent the positions of the atomic layers at the interfaces. 
The position z is measured from the top metal atom. 



0.5 



o 



^ 



-0.5 



.v.*;- 

£' . S 


1 ^ 
/ \ 

h \ \ 

w 1 \ 




- V/A 

1/16 V/A 
2/16 V/A 
3/16 V/A 
4/16 V/A 

- 5/16 V/A 

- 6/16 V/A 

- 7/16 V/A 










t ' ■ ' • '\ 

t • ."A 
/ S 


HtT — \~\ 














: \ '"' 
\ \\' 1 




•>7-BN 


V-' G 


\' 




\ 


■ -1+ 







5 

z(A) 



10 



FIG. 9. (color online) Influence of an external electric field 
on the BN|Gr interface dipole layer. The plane-averaged elec- 
tron density difference AnBN|Gr('2;5 ^ext) is shown for different 
electric field strengths (solid lines). The dotted and dashed 
lines show the polarization of isolated /i-BN and Gr sheets, 
respectively, in these external fields. The vertical dashed lines 
represent the positions of the atomic layers at the interfaces. 



Ancu|BN(^7 ^ext) for /i-BN adsorbed on the Cu(lll) sur- 
face calculated with £;ext = -0.44, and +0.44 V/A. 
This fi gure clearly demonstrates that A77/Cj|bn(z, .£/ext) 
is essentially independent of the electric field strength 
within this range, and therefore also the interface dipole 
and the potential step. A similar independence can 
be observed in the electronic displacement at a h- 
BN|graphene interface in Fig. 9 (solid lines). Hence 
the interface dipoles, and the potential steps Am|bn and 
^BN|Gr are truly intrinsic properties of these interfaces. 

The fact that the electronic displacement at the in- 
terfaces is field-independent (for the field strengths used 
here) does not mean that the electron densities are so. 




3 6 

z(A) 



FIG. 8. (color online) Influence of an external electric field 
on the Cu(lll)|BN interface dipole layer. The plane- averaged 
electron density difference Ancu|BN('2;5 ^ext) is shown for three 
different electric fields (perpendicular to the interface). The 
vertical dashed lines represent the positions of the atomic lay- 
ers at the interface. 



as both /i-BN and graphene are polarized in an elec- 
tric field. Fig. 9 shows the changes in electron density 
n(z,£^ext) — ^(^,0) upon applying an electric field, both 
for an isolated graphene layer (dashed lines) and for an 
isolated /i-BN slab (dotted lines). The changes clearly 
illustrate the induced dipoles on the individual layers. 

For a complete met al|/i-BN| graphene stack, we visual- 
ize the total electronic displacement caused by interface 
formation, as well as by the external electric field, as 

Antot(^,£^ext) = 

^M|BN|Gr(^, ^ext) " ^m(^) - ^Bn(^) - ^Gr(^), (15) 

where the reference charge densities nM, ^bn, ^Cr are for 
the isolated systems in zero external field. Results are 
shown in Fig. 10 for a Cu(lll)|/i-BN|graphene structure 
with five layers of /i-BN. The dipoles at the Cu(lll)|/i- 
BN and /i-BN| graphene interfaces are essentially the 
same as the dipoles at individual isolated interfaces, see 
Figs. 8 and 9, implying that the two interfaces in the 
Cu(lll)|/i-BN|Gr stack are not coupled. 

Antot in Fig. 10 for the zero external field case (£/gxt — 
0) shows small oscillations in the /i-BN slab, indicat- 
ing that the /i-BN layers are polarized. The effective 
work function difference between the Cu substrate and 
graphene, Eq. (8), results in electron transfer from Cu to 
graphene across the /i-BN slab to establish a single Fermi 
level in equilibrium. This sets up an intrinsic electric 
field in the /i-BN slab, see Eq. (7). This charge transfer 
can be counteracted by applying an external electric field 
£;ext = -£^0, see Eq. (6). 

This external field results in charge neutral (i.e. un- 
doped) graphene, but the field in the /i-BN slab is still 
non-zero, see Eq. (4) with a = and £^ext = — ^o- An ex- 
ternal field of £;ext = £^io ^ -0.43 V/A, Eq. (9), gives a 
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FIG. 10. (color online) Plane- averaged electron density differ- 
ence Antot(^, ^ext) for a Cu(lll)|5-layers-of-/z-BN|graphene 
structure for external electric fields Ee^t — — ^o,0,^io V/A. 
The vertical dashed lines represent the positions of the atomic 
layers at the interfaces. 



zero internal field, see Fig. 10. Note that, whereas Antot 
in the region between the topmost /i-BN layer and the 
graphene sheet is almost field-independent, it strongly 
depends on the field on the graphene side {z > 30 A). 
This just reflects the polarization behavior of the isolated 
graphene sheet, see Fig. 9. 



Fermi level shifts in met al|/i-BN| graphene 
heterostructures 



All the parameters that appear in Eqs. (6-8) have now 
been calculated by studying constituent parts of M|/i- 
BN|Gr stacks. Table I, and the model can be tested by 
comparing its predictions to the results obtained from 
explicit DFT calculations on full M|/i-BN|Gr structures. 
Fig. 6 demonstrates that the electric field dependence of 
the Fermi level shift is described very well by the model. 
We would like to emphasize that these curves do not 
represent data fits, as all parameters in the model are 
fixed. 

We next study how the thickness of /i-BN, measured 
in terms of the number of discrete layers n, influences the 
Fermi level shift A^p in graphene. In Table II, AE^ cal- 
culated with DFT for Cu( 1 11) |nx/i-BN| graphene struc- 
tures with < n < 6 layers is compared to values 
predicted by the model for a vanishing external field, 
A^^. The agreement between the two sets of results is 
to within 12 meV or better. The n = case corresponds 
to graphene directly adsorbed onto the Cu(lll) surface 
and the results agree with those found previously.^'^ The 
good agreement found in Table II indicates that the two 
interfaces are not strongly coupled. Even for a single 
/i-BN layer, n = 1, the model works well. 

In the absence of an applied potential, graphene in a 
Cu(lll)|/i-BN|graphene stack is doped n-type in equi- 



TABLE II. Calculated work function W and Fermi energy 
shift A^F of Cu(lll)|nx/i-BN|graphene structures with n 
layers of /i-BN, and Fermi energy shifts A£^™ predicted by 
the model. All values are in eV. 



n 





1 


2 


3 


4 


5 


6 


W 

AEf 

AE^ 


4.386 
-0.171 
-0.17^ 


4.223 
-0.304 
-0.297 


4.299 
-0.234 
-0.246 


4.325 
-0.206 
-0.216 


4.336 
-0.193 
-0.196 


4.346 
-0.180 
-0.181 


4.354 
-0.171 
-0.169 



"Ref. 5 



librium. The degree of doping decreases as the thick- 
ness of the dielectric layer is increased. This can be 
simply understood from Eq. (6) where for large values 
of the Fermi level shift scales asymptotically with d 
as AEf ^ 1/V^. Fig. 11 illustrates this scaling, and 
also demonstrates the deviation of this scaling for small 
d. The good agreement between the results of the first- 
principles calculations and the model seems to improve 
with increasing /i-BN thickness. 

As discussed in the previous section, the Fermi level 
can be positioned at the charge neutrality point of 
graphene by application of an external electric field 
^ext = —Eq. According to Eq. (7), Eq should scale with 
the /i-BN layer thickness d diS Eq (x 1/d. This scaling is 
demonstrated in Fig. 12. The symbols represent the £^ext 
values needed in the DFT calculations to make graphene 
charge neutral in stacks with variable /i-BN thickness. 
From Eq. (11) we then predict that the gate voltage re- 
quired to shift the Fermi level to the graphene conical 
points should be independent of how thick /i-BN is. 

The dependence of the doping level, both as a function 
of the external field £^ext7 and as a function of the h- 
BN layer thickness, is shown in Fig. 13. The agreement 
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FIG. 11. (color online) The Fermi level shift versus the thick- 
ness d of the /i-BN film, plotted as a function of l/\/d and 
(inset) as a function of d. The symbols are the results of DFT 
calculations on Gu(lll)|/i-BN| graphene structures with 1-6 
layers of /i-BN. The dashed line represents the model calcu- 
lated with Eq. (6), the dotted line represents the asymptotic 
l/Vd dependence. 
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FIG. 12. (color online) The value of the intrinsic electric field 
Eq (Eq. (7)) as a function of 1/d where d is the thickness of the 
/i-BN film. An external electric field E^ext — —Eq is required 
to restore graphene to charge neutrality. The symbols are 
the results of DFT calculations for Cu(lll)|/i-BN|graphene 
structures with 1-6 layers of /i-BN. 



0.4 



0.2 



> 



-0.2 



-0.4 




Cq.^^ Ag 

-0.5 0.5 

(eV) 



-1 



FIG. 14. (color online) The Fermi level shift for metal 
BN I graphene structures with 2 layers of /i-BN, plotted against 
Vb, according to Eqs. (12) and (8) with = and parameters 
from Table I. The line is for a fixed d, and the circles are for 
the equilibrium metal|/i-BN bonding distances. 



between the model and the DFT results is very good. 
The characteristic shape of the curves is determined by 
Eq. (6). The curves for different /i-BN layer thicknesses 
cross at £^ext ^ —0.4 V/A. At this external field strength 
the doping level is independent of the thickness of the h- 
BN dielectric spacer. This point corresponds to the field 
^io = —0.43 V/A discussed above, where the internal 
field in the /i-BN spacer is zero. 

For a dielectric consisting of two /i-BN monolayers, 
Figure 14 shows how graphene in a metal |/i-BN| graphene 
structure is doped in the absence of an external field by 
different metals. Different metals have different work- 
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FIG. 13. (color online) The Fermi level shift versus external 
electric field strength for a Cu(lll)|/i-BN| graphene structure 
with 2-5 layers of /i-BN. The lines represent the model of 
Eq.(6) for different thicknesses of dielectric. The symbols are 
values of the doping level calculated with DFT. 



functions Wu and form different interface dipoles with 
the dielectric; each metal can be characterized by a dif- 
ferent value of Vq given by Eq. (8) using the parameters 
given in Table I. The dashed line represents the model 
of Eq. (12) for a fixed distance the circles correspond 
to the different metals where the equilibrium distances of 
Table I have been used. The predicted Fermi level shift 
AEp varies with the metal used for the gate electrode 
from +0.18 eV for Pt to -0.39 eV for Co. 

Most of the metals considered give rise to n-type dop- 
ing of graphene; p-type doping is only found for the high 
work function metals Au and Pt. The amount of doping 
is however much less than what would be expected on 
the basis of the work function difference between these 
metals and graphene only. This is because of the large 
potential step Am|bn between Au or Pt and /i-BN, see 
Table I, which effectively lowers the metal work function. 
Although most of the other metals have a work function 
that is comparable to or higher than that of graphene, 
they result in n-type and not in p-type doping. Again 
this can be attributed to the effect of the potential step 
Am|bn- a similar effect was found for graphene that is 
directly adsorbed on a metal substrate. ^'^ 

An external electric field can be used to counter- 
act the intrinsic doping of graphene in these metal|/i- 
BN I graphene stacks and position the Fermi level at the 
charge neutrality point of graphene corresponding to 
AEp = 0. The field strengths required for the differ- 
ent metals cover a range of ~ 0.5 V/A. For the stacks 
considered in Fig. 14, this corresponds to a range of 2 V 
in the gate voltage. 
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V. SUMMARY AND CONCLUSIONS 

We have developed a model, Eq. (6), for the electro- 
static doping of graphene in a field-effect device that 
describes the position of the Fermi level in graphene 
for a metal I dielectric I graphene structure as a function of 
the gate voltage. As input parameters the model needs 
the work function of the metal, the thickness and the 
dielectric constant of the dielectric material, and the 
potential steps formed at the metal | dielectric and the 
dielectric I graphene interfaces. All of these parameters 
can be obtained from DFT calculations on single surfaces 
or interfaces, which makes the model parameter- free. 

The model predicts doping of graphene even when the 
bias voltage is zero. This effect is caused by charge equi- 
libration, i.e. the need to establish a common Fermi level 
across the metal | dielectric | graphene stack. The amount 
of doping is not only guided by the difference in work 
function between the metal and graphene, but also to 
a significant extent by the potential steps at the inter- 
faces. These potential steps result from interface dipole 
layers that are formed by the chemical interactions at 
the interface. In particular the potential step at the 
metal I dielectric interface can be large, i.e., > 1 eV, even 
when the bonding between the metal and the dielectric 
material is not strong. This potential step effectively low- 
ers the work function of the metal, which causes metal 
substrates with work functions < 5 eV to dope graphene 
n-type. Only higher work function metal substrates lead 
to p-type doping. 

The amount of doping at zero bias decreases with 
increasing thickness of the dielectric. The Fermi level 
shift in graphene as a function of the applied electric 
field resulting from the gate voltage has a square-root 
like behavior. A similar behavior has been observed 
in experiments on gated p-Si|Si02| graphene structures, 
by work function measurements,^^ and by conductivity 
measurements.^^ According to our model, this square- 
root like behavior is a direct result of the density of states 
in graphene being linear around the conical points and 
implies that the adsorption of graphene on the dielectric 
in those experiments was sufficiently weak as to leave the 
electronic structure of graphene essentially unchanged. 

The model is tested by comparing its predictions 
to results obtained from explicit DFT calculations on 
metal I dielectric I graphene structures in the presence of an 
external electric field. We use the (111) or (0001) sur- 
faces of the fee or hep close-packed metals Al, Co, Ni, 
Cu, Pd, Ag, Pt, Au, which span a range of work func- 
tions from 4.2 to 6.0 eV. As our dielectric we use either 
vacuum or /i-BN, which is an insulating material with a 
similar honeycomb structure as graphene, and we vary 
the thickness of the dielectric by varying the number of 
/i-BN layers. The model gives a good description of the 
Fermi level shift in graphene as a function of the applied 
electric field even for a /i-BN dielectric only a monolayer 
thick. The thickness dependence of the Fermi level shift 
is also reproduced very well by the model. 



To be able to perform microscopic DFT calculations 
it was necessary to make approximations which are in 
general not satisfied by real materials. In particular, we 
assumed that graphene is commensurate with /i-BN and 
that both are commensurate with the underlying metal 
substrate. Where the in-plane bonding is very strong 
and the inter-plane bonding is weak as is the case for 
graphene on /i-BN, the resulting composite system re- 
mains incommensurate^^'^^'^^'^^ and it is not possible to 
carry out first-principles electronic structure calculations. 
Even when the bonding of graphene or /i-BN to a metal 
substrate is stronger, the periods that are found for com- 
mensurability are frequently so large^^ that DFT calcu- 
lations are exorbitantly expensive. In such cases, we sug- 
gest that the parameters in the model we have outlined 
should be regarded as free parameters and determined 
by fitting to experimental observations. In this way, it 
should be possible to apply the model to very much more 
complex systems such as Si02 or SiC dielectrics and poly- 
crystalline Si gates. Intrinsically planar dielectric spacers 
like M0S2 and WS2 that are currently of interest because 
of their potential in graphene electronics^^'^^ are inter- 
mediate; they can be studied using DFT calculations, 
modelling the interfaces with metals and graphene using 
large lateral supercells, or the relevant parameters can be 
determined by fitting to experiment, or a combination of 
both. 
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Appendix: Graphene Brillouin zone sampling 

The first-principles calculations presented in this paper 
were carried out with a 36 x 36 sampling of the two- 
dimensional (2D) Brillouin zone (BZ) for (1 x 1) graphene 
or an equivalent sampling density for (2x2) graphene. At 
energies close to the charge neutrality level of graphene, 
the density of states (DoS) is described well by the linear 
expression D{E) = Do\E\/A with A the surface area of 
a graphene unit cell. Table III shows how the parameter 
1^0 depends on the number of k-points used to sample the 
2D-BZ. An unexpectedly dense grid (288 x 288 points) is 
required to obtain a converged DoS which in view of the 
rather simple linear behavior of the DoS is at first glance 
surprising. 

The main reason for requiring such a dense grid lies in 
the fact that for energies within 1 eV of the Fermi level. 
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Do states/ (eV^ unit cell) 





Freestanding 


Graphene 


Graphene 


on h-BN 


N 


^ < 


^ > 


E <0 


E>0 


36 


0.088 


0.092 


0.102 


0.102 


72 


0.112 


0.119 


0.121 


0.112 


144 


0.109 


0.116 


0.114 


0.115 


288 


0.108 


0.115 


0.114 


0.114 



TABLE III. Values of Do calculated using a grid obtained 
by dividing the reciprocal lattice vectors of the graphene lat- 
tice into N intervals, and using the tetrahedron method^^ to 
calculate the DoS. All grids contain the K-point. 



all of the states contributing to the DoS lie in a very 
small part of the BZ around the conical points. This is 
illustrated in Fig. 15 where it can be seen that even a 
36 X 36 k-point grid leads to a very coarse sampling of 
the relevant bands close to the conical points. 

The DoS for various grid densities is shown in Figure 
16, calculated with the linear tetrahedron method, as 
implemented in VASP.^^^'^ For the 36 x 36 grid the DoS 
is perfectly linear up to ±0.45 eV, where it jumps to a 
different value. These jumps occur at eigenvalues cor- 
responding to the k-points that are explicitly included 
in the grid. With increasing grid density the number of 
jumps increase, but their size decreases, such that the 
DoS converges to a linear function. 

The jumps in the DoS at the k-grid points are not 
an effect of the trigonal distortion of the graphene Dirac 
cone; the DoS of an ideal conical dispersion E (x \k\ (cen- 
tered on r) shows similar jumps when calculated with the 
tetrahedron method. The slopes I^o listed in Table III 
were extracted by fitting a linear curve to the DoS for 
energies +0.4 (or —0.4) eV with respect to the charge 
neutrality energy. For free-standing graphene the con- 
verged LDA values are 0.108 and 0.115 /(eV^ unit cell) 
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FIG. 16. (color online) DoS of graphene close to the Dirac 
point obtained by successively doubling the linear sampling 
density of the k-point grid used in the tetrahedron scheme. 
With increasing grid density the DoS converges to a smooth 
line. 



for the filled and empty states, respectively. 

The graphene DoS also depends somewhat on the sub- 
strate used, as illustrated by Fig. 17 and Table III. The 
slope 1^0 for graphene adsorbed on /i-BN is somewhat 
larger than that of free-standing graphene for holes £^ < 
and somewhat smaller for electrons E < ^. Curiously, 
the DoS is more symmetric for electrons and holes for 
graphene on /i-BN. A larger means that the graphene 
TT-bands, and therefore the Dirac cone, become somewhat 
flattened, as illustrated in Fig. 18. 

The DFT calculations on M|BN|Gr structures reported 
in this manuscript were carried out using a 36 x 36 k-point 
grid. For internal consistency, a value I^o = 0.102/(eV^ 
unit cell) corresponding to this grid (see Table III) was 
used in the model that fits the numerical results essen- 




FIG. 15. (color online) a) Irreducible Brillouin zone (IBZ) of 
the graphene unit cell sampled with a 6 x 6 k-point grid, b) 
The eigenvalues of free-standing graphene in the vicinity of 
the Dirac (K) point sampled with a 36 x 36 k-point grid. The 
dashed contour about the K point in a) and b) represent a 
cut through the Dirac "cone" at -1.0 eV from the Dirac point 
while the solid contour in b) represents a cut at -0.5 eV. 




E-Ep (eV) 

FIG. 17. (color online) DoS of free-standing graphene 
(red/gray) and of graphene on /i-BN (black), close to the 
Dirac point. The solid and dashed curves were calculated 
using 36 X 36 and 288 x 288 grids, respectively. 




FIG. 18. (color online) Cross sections of the graphene Dirac 
cone at energies —0.1, —0.5, —1.0 eV with respect to the 
charge neutrality point, as calculated using a 288x288 k-point 
grid. The black lines correspond to graphene on /i-BN and 
the red (gray) lines to free-standing graphene. The axes are 
in units of ^tt/^. 
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